The virial theorem and the dark matter problem in hybrid metric-Palatini gravity 
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Hybrid metric-Palatini gravity is a recently proposed theory, consisting of the superposition of the 
metric Einstein-Hilbert Lagrangian with an f(7Z) term constructed a la Palatini. The theory predicts 
the existence of a long-range scalar field, which passes the Solar System observational constraints, 
even if the scalar field is very light, and modifies the cosmological and galactic dynamics. Thus, the 
theory opens new possibilities to approach, in the same theoretical framework, the problems of both 
dark energy and dark matter. In this work, we consider the generalized virial theorem in the scalar- 
tensor representation of the hybrid metric-Palatini gravity. More specifically, taking into account the 
relativistic collisionless Boltzmann equation, we show that the supplementary geometric terms in the 
gravitational field equations provide an effective contribution to the gravitational potential energy. 
We show that the total virial mass is proportional to the effective mass associated with the new terms 
generated by the effective scalar field, and the baryonic mass. This shows that the geometric origin 
in the generalized virial theorem may account for the well-known virial theorem mass discrepancy 
in clusters of galaxies. In addition to this, we also consider astrophysical applications of the model 
and show that the model predicts that the mass associated to the scalar field and its effects extend 
beyond the virial radius of the clusters of galaxies. In the context of the galaxy cluster velocity 
dispersion profiles predicted by the hybrid metric-Palatini model, the generalized virial theorem can 
be an efficient tool in observationally testing the viability of this class of generalized gravity models. 
Thus, hybrid metric-Palatini gravity provides an effective alternative to the dark matter paradigm 
of present day cosmology and astrophysics. 

PACS numbers: 04.50.+h, 04.20.Jb, 04.20.Cv, 95.35.+d 



I. INTRODUCTION 

Modern astrophysics and cosmology are facing two ma- 
jor intriguing challenges, namely, the dark energy and 
the dark matter problems, respectively. Two important 
astronomical observations, the fiat rotations curves of 
galaxies, and the virial mass discrepancy in clusters of 
galaxies led to the necessity of inferring the existence 
of a special form of matter, called dark matter, which 
interacts with baryonic matter only gravitationally, and 
whose presence can explain the dynamical behavior of 
test particles at the galactic and extra-galactic scale. On 
the galactic or intra-galactic scale, the astronomical ob- 
servations show a linear mass increase even in regions 
where very little baryonic matter can be detected, with 
the rotational velocities attaining an approximately con- 
stant value, Vtgoa ~ 200 — 300km/s, within a distance r 



'Electronic address: capozzie@na.infn.it 
^Electronic address: tiberiu.harko@gmail.com 
tElectronic address: tomi.koivisto@fys.uio.no 
^Electronic address: flobo@cii.fc.ul.pt 
^Electronic address: gonzalo.olmo@csic.es 



from the center of the galaxy [l| . 

The determination of the total mass that is gravita- 
tionally bounded in clusters of galaxies is important for 
our understanding of the nature and evolution of struc- 
tures on cosmological scales. The gravitational masses of 
clusters of galaxies are estimated by assuming a hydro- 
static equilibrium of both the hot intra-cluster gas and of 
the galaxies with the binding cluster potential. Therefore 
the total mass of a cluster of galaxies can be estimated 
in two ways. In the first method, by taking into account 
the dynamical motions of the member galaxies of the 
cluster, and with the application of the virial theorem, 
one obtains an estimate My for the mass of the cluster. 
Secondly, the total baryonic mass Mb can be determined 
by adding the mass of each individual galaxy member of 
the cluster. The mass discrepancy at the galactic cluster 
level arises as observations show that My is much greater 
than M B , with typical values of My jM - 20 - 30 Q. 

These important astrophysical observations are usu- 
ally explained by postulating the existence of a new form 
of matter, called dark matter, assumed to be a cold 
pressure-less fluid extended in a spherically symmetric 
halo around the galaxies and clusters (for a recent re- 
view on the dark matter properties see Q). Many candi- 
dates for dark matter particles have been proposed, the 
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most popular ones being the Weakly Interacting Massive 
Particles (WIMPs) !3(. The interaction cross section of 
WIMPs with standard baryonic matter is assumed to be 
extremely small, but non-zero, and therefore these parti- 
cles may be directly observable experimentally. Recently, 
a pressure-less and non-comoving two-fluid dark mat- 
ter model has also been analysed, in which dark matter 
is represented as a two-component fluid thermodynamic 
system, without interaction between the constituent par- 
ticles of different species [1] . 

However, despite the extensive experimental and ob- 
servational investigations, up to now non- gravitational 
evidence for dark matter is still lacking. Moreover, ma- 
jor accelerator and reactor experiments, like LHC, did 
not bring convincing evidence for the existence of new 
physics beyond the standard model, on which the dark 
matter hypothesis is based upon. Therefore, the possi- 
bility that Einstein's (and Newton's) theory cannot de- 
scribe gravitational phenomena at the scale of galaxies 
and clusters of galaxies must not be excluded a priori. 
Several theoretical approaches, in which "dark matter" 
can be understood as a modification of the gravitational 
laws at large scales have been extensively proposed in the 
literature [5|. 

The virial theorem, which gives a simple relation be- 
tween the kinetic and potential energy of a system of 
particles, plays an important role in astrophysics and cos- 
mology 6]. By observing the velocities of test particles 
and by assuming hydrostatic equilibrium, with the use of 
the virial theorem, one can obtain the mean density of 
astrophysical objects such as galaxies, clusters, or super 
clusters. An important application of the virial theorem 
is the determination of the total mass of the clusters of 
galaxies. The virial theorem is also a powerful tool for 
the study of the stability of gravitationally bounded ob- 
jects. The extension of the Newtonian virial theorem to 
the general relativistic case has led to several versions of 
the virial theorem m , including the effect of a cosmo- 
logical constant [1, [j| , the generalization to brane world 
models [T(J, f(R) gravity jnlLDGP brane models [12], 
and to Palatini f(R) models [13j . 

In this context, a novel approach to modified theories 
of gravity was proposed [lj], that consists of adding to 
the Einstein-Hilbcrt Lagrangian a f(R) term constructed 
within the framework of the Palatini formalism. Using 
the respective dynamically equivalent scalar-tensor rep- 
resentation, even if the scalar field is very light, the the- 
ory passes the Solar System observational constraints. 
Therefore the long-range scalar field is able to modify the 
cosmological and galactic dynamics, but leaves the Solar 
System unaffected. The absence of instabilities in per- 
turbations was also verified, and explicit models, which 
are consistent with local tests and lead to the late-time 
cosmic acceleration were also found. 

The cosmological applications of the hybrid metric- 
Palatini gravitational theory were investigated in (l5j . 
Criteria to obtain cosmic acceleration were analyzed, and 
the field equations were formulated as a dynamical sys- 



tem. Several classes of cosmological solutions, depending 
on the functional form of the effective scalar field poten- 
tial, describing both accelerating and decelerating uni- 
verses, were explicitly obtained. Furthermore, the cosmo- 
logical perturbation equations were derived and applied 
to uncover the nature of the propagating scalar degree 
of freedom and the signatures of these models predicted 
in the large-scale structure. In addition to this, the hy- 
brid metric-Palatini theory was considered in wormhole 
physics [l6j]. The general conditions for wormhole solu- 
tions according to the null energy conditions at the throat 
in the f(X) gravity were presented in [l6j . Several worm- 
hole type solutions were also obtained and analyzed. In 
the first solution, the redshift function and the scalar 
field were specified, and the potential was chosen so that 
the modified Klein-Gordon equation can be simplified. 
This solution is not asymptotically flat and needs to be 
matched to a vacuum solution. In the second example, 
by adequately specifying the metric functions and choos- 
ing the scalar field, one can find an asymptotically flat 
spacetime. 

The purpose of the present paper is to consider the 
virial theorem in the framework of hybrid metric-Palatini 
gravity. Using the collisionless Boltzmann equation in 
the modified Einstein field equations we derive a gener- 
alized virial equality for the hybrid model. The gener- 
alized virial theorem takes into account the presence of 
the supplementary scalar field related terms, which do 
appear due to the modification of the gravitational ac- 
tion. These supplementary geometric terms, and their 
scalar field equivalent, give an effective contribution to 
the gravitational potential energy, with the total virial 
mass being given by the sum of the effective mass as- 
sociated to the new scalar field related terms, and the 
baryonic mass, respectively. Therefore the new scalar 
field may account for the virial theorem mass discrepancy 
in clusters of galaxies The gravitational field equa- 
tions of the hybrid metric-Palatini gravitational model 
together with the virial theorem also allow to obtain the 
metric inside the cluster of galaxies in a simple form, as 
functions of physical parameters that can be fully deter- 
mined from astrophysical observations, like, for example, 
the temperature of the intra-cluster gas and the radius 
and central density of the cluster core. Therefore the gen- 
eralized virial theorem in hybrid metric-Palatini gravity 
can be an efficient tool in observationally testing the vi- 
ability of this class of generalized gravity models. 

The present paper is organized as follows. In Sec- 
tion HH we briefly present the general formalism of hy- 
brid metric-Palatini gravity, for self-completeness and 
self-consistency. In Section IIII1 the gravitational field 
equations for spherically symmetric galactic clusters in 
the scalar-tensor version of hybrid metric-Palatini grav- 
ity and the relativistic Boltzmann equation are explored. 
In Section IIV1 the generalized virial theorem in hybrid 
metric-Palatini gravity is derived. Astrophysical appli- 
cations of the virial theorem are explored in Section [V] 
In particular, predictions of the geometric mass and geo- 
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metric radius from galactic cluster observations are pre- 
sented, and the behavior of the galaxy cluster velocity 
dispersion in hybrid metric-Palatini gravity models is 
considered. Finally, in Section IVI1 we discuss and con- 
clude our results. 



II. HYBRID METRIC-PALATINI GRAVITY: 
THE FORMALISM 

The action for hybrid metric-Palatini gravity is given 
by 



where C m is the matter Lagrangian density. The scalar 
curvature ft depends on both the metric and an indepen- 
dent dynamical connection T" v , and is defined as 



(2) 

where 1Z^ V is the Ricci tensor obtained from the connec- 
tion f^. 

Varying the action (|T|) with respect to the metric, one 
obtains the following gravitational field equations 

+ F(K)K^ - ^f(K)g^ = k 2 T^ , (3) 



where the matter stress-energy tensor is defined as 
2 S(y^gC m ) 



f=g Sg^ 



(4) 



The independent connection is compatible with the 
metric F{TZ)g^ u , conformal to g^ u ; the conformal factor 
is given by F(1Z) — df(lZ)/dlZ. The latter considerations 
imply that 



ft/^ RflU 



1 



" v ' 2F 2 (ft) 
1 



F(n) 



F{n) tfl F{Tl), v 
1 1 



2F(ft) 



g^DFCR.) . (5) 



Note that ft can be obtained from the trace of the field 
equations ©, which yields F(ft)ft - 2/ (ft) - R = k 2 T. 

Introducing an auxiliary field, the hybrid metric- 
Palatini action ([T]) can be turned into a scalar-tensor the- 
ory given by the following action (we refer the reader to 
for more details) 



s = i / dix ^9 [R + <pn- v{cj>)] + s„ 



(6) 



Varying this action with respect to the metric, the scalar 
cj> and the connection yields the following field equations 

R^ + ^Tl^- -(R + <pn-V) gtlu = k 2 T^, (7) 



ft- V 



o, 

0, 



(8) 
(9) 



respectively. Note that the solution of Eq. ([9]) implies 
that the independent connection is the Levi-Civita con- 
nection of a metric h^ v = <f>g^ u - Thus we are dealing 
with a bi-mctric theory and ft M! , and R^ v are related by 
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ft M y = R^v + 202 ^M^ 7 *- 



- I V M V„ V 



(10) 

and consequently ft = R + 3V A1 0V /i < ?i/2(/) 2 - 
which can be used in the action © to get rid of the 
independent connection and obtain the following scalar- 
tensor representation: 



S=^- 2 fd^^-9[R + f(n) + 2n 2 C m ] , (1) s= 1 f di 
J 2k 2 J 



(11) 

where S m is the matter action. Note that the hy- 
brid metric-Palatini gravitational theory is equivalent 
with the purely metric Brans-Dicke-like action, with the 
Brans-Dicke parameter w — —3/2, but with a different 
coupling to matter. 

The variation of this action with respect to the metric 
tensor gives the field equations 



(l + flGfu, = k 2 T m „ + V„V„<; 



_3_ 

U 



1^ - -Vg^ v ,{l2) 



Varying the action with respect to the scalar field we 
obtain 
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24>- 



dV 







Moreover, one can show that the identity 

,dV 



2V 



k 2 T + R 



(13) 



(14) 



also holds, and that the scalar field 4> is governed by the 
second-order evolution equation 



□« 



1 

24 



1 



IV - (1 



dV 



which is an effective Klein-Gordon equation. 



(15) 



III. GALACTIC CLUSTERS: FIELD 
EQUATIONS FOR A SYSTEM OF IDENTICAL 
AND COLLISIONLESS POINT PARTICLES 

Astronomical and astrophysical observations have 
proved that galaxies tend to concentrate in larger struc- 
tures, called clusters of galaxies, with total masses rang- 
ing from 10 13 M Q for groups up to a few 10 15 M Q for 
very large systems. The cluster morphology is usually 
dominated by a regular, centrally peaked main compo- 
nent [HI, [l9| . Since clusters are considered to be "dark 
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matter" dominated astronomical systems, their forma- 
tion and evolution is controlled by the gravitational force. 
The mass function of the clusters is determined by the 
initial conditions of the mass distribution set in the early 
universe. The evolution of the cosmological perturba- 
tions of the large scale matter distribution on scales com- 
parable to the size of the clusters is linear. The gravi- 
tational growth of the cosmological density fluctuations 
and the gravitational instabilities leading to cluster for- 
mation have been extensively studied by using analytical 
methods, as well as numerical simulations |20|. 

In order to derive the generalization of the relativistic 
virial theorem for galaxy clusters in the hybrid metric- 
Palatini gravitational models we need, as a first step, to 
obtain the gravitational field equations for a static and 
spherically symmetric distribution of matter. To thus 
effect, consider a self-gravitating system of identical, col- 
lisionless point particles in random motion. To obtain the 
basic field equations we will use the scalar-tensor repre- 
sentation of hybrid metric-Palatini gravity, which allows 
a clear physical interpretation of the model. 

We assume that the geometry of the cluster can be de- 
scribed by a time-oriented Lorentzian four-dimensional 
space-time manifold M. . The metric of an isolated spher- 
ically symmetric cluster is given by 



ds 2 



Mr) dt 2 + e X(r) dr 2 + , 



+ sin 2 9d<p 2 ) . (16) 



The galaxies in the cluster are considered identical and 
collisionless point particles, and their space-time distri- 
bution is described by a distribution function fg. The 
latter function obeys the general relativistic Boltzmann 
equation, which will be presented in the next Section. 

Thus, the energy-momentum tensor of matter is de- 
termined by the distribution function /g, and is given 

by [13 



/ 



Tuu = / Jbtti u^Uv du, 



(17) 



where m is the mass of the particle (galaxy), m m is the 
four-velocity of the galaxy and du = durdugdu^/ut is 
the invariant volume element of the velocity space. The 
energy-momentum tensor T^ v of the matter in a cluster 
of galaxies can be represented in terms of an effective 
density p e s and of two effective anisotropic pressures, the 
radial and the tangential p^g pressures, respectively, 
given by 

Pcff = p(u 2 t ) ,pXi = p(u 2 r ) ,p[^ = p(u 2 e ) = p{u%) , 

(18) 

where p is the mass density of the ordinary baryonic mat- 
ter, and (u 2 ), i = t,r,d,<p is the average value of uf, 
i = t,r, 0jjp, the square of the components of the four- 
velocity d. 

By using this form of the energy-momentum tensor, 
the gravitational field equations describing a cluster of 



galaxies in hybrid metric-Palatini gravity take the form 



- e -J + ~o = -p (u 2 ) 

\r 2 r J r 2 l + 4> x *' 

' -V^ + ^—^^-U) ' 



1 



2(1 

' ' -V^VV+t ±l }, ^ V X 6V X 6, (19) 



26(1 



46(1 + 6) 



imT6) V ^ ~ 7TmV*V**V^ (20) 



46(1 + 6) 



G 



-t)=^t+-A<) 



+ WW) vi<p) -—6^ ve -^ 



3 1 



2 6(1 



-Vg6V b 



3 1 



46(1 



-V X 6V X 6(21) 
G 



- + 2(TT0)^ W TT0 (V ^- D) ' 
i — L_ v ^v^ 3 1 



2 6(1 + 6) 



46(1 + , 



V A ^V A ^22) 



where there is no summation upon the pair of indices 
(t, r, 9, ip). A useful relationship is obtained by adding the 
gravitational field equations Eqs. ([T9"l) -(f2"2" |) . from which 
we obtain the following equation 



—A 



v" v 12 v' v'\' 



G 



(23) 



where (u 2 ) = (u 2 ) + (u 2 ) + (u 2 ) + (u%). 

Since we are interested in astrophysical applications at 
the extra-galactic level, we may assume that the devia- 
tions from standard general relativity (corresponding to 
the background value 6 — 0) are small, i.e., 6 <C 1. Thus, 
Eq. (|23|) can be rewritten as 



v" v 12 v> v'\ . , 
< ^ + — + -- — \- ^Gp(u 2 ) + 4^G P: 



(eff) 

(24) 



where 



4vG P f n ~V(6) + (2V t V* +U)6- -V t 6V l 6, (25) 

represents an effective energy of the scalar field in the 
hybrid metric-Palatini gravitational model. 
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IV. THE GENERALIZED VIRIAL THEOREM 
IN HYBRID METRIC-PALATINI GRAVITY 

As a second step in the derivation of the virial theorem 
for galaxy clusters, which we describe by using the distri- 
bution function f B , we have to write down the Boltzmann 
type equation governing the evolution of the distribution 
function. This equation can then be integrated over the 
velocity space, to yield an equation which, used together 
with the gravitational field equations, gives finally the re- 
quired generalization of the virial theorem for a general 
relativistic distribution of point particles. 



A. The relativistic Boltzmann equation 

The transport equation for the distribution function 
for a system of particles in a curved arbitrary Riemannian 
space-time is given by the Boltzmann equation 17 1 



^-pVit.A ]/, = ... 



(26) 



where p a is the four-momentum of the particle, and T l a p 
are the Christoffel symbols associated to the metric. Note 
that the collissionless Boltzmann equation states that 
the local phase space density viewed by an observer co- 
moving with a star or galaxy is conserved. 

A considerable simplification of the mathematical for- 
malism can be obtained by introducing an appropriate 
orthonormal frame or tetrad e^(x), a = 0,1,2,3, which 
varies smoothly over some coordinates neighborhood U 
and satisfies the condition /"e'ej = r/' 
where i] ab is the Minkowski metric tensor 
tangent vector p^ at x can be expressed as p^ — p"e 
which defines the tetrad components p a . 

In the case of the spherically symmetric line element 
given by Eq. (fl~6|) an appropriate choice of the frame of 
orthonormal vectors is [8|, [l7| : 



6 for all x € U, 
' ll- Any 
a ■ 



= e^ 2 S° 



e\ = e x ' 2 5 



e 2 = rS z 



ef, = r sin 98f, . 



(27) 

Let u M be the four-velocity of a typical galaxy, satis- 
fying the condition u^u^ = —1, with tetrad components 
u a = u^e^. In tetrad components the relativistic Boltz- 
mann equation is 



u a e" 



df B 



■ lbc U U 



, df B 
du a 



= 0, 



(28) 



where the distribution function f B = fs{x ll ,u a ) and 
Ibc = e ti;v e b e c are tne Ricci rotation coefficients \§L \lT\. 
By assuming that the only coordinate dependence of the 
distribution function is upon the radial coordinate r, 



Eq. ([28)1 becomes 

df B 
ui-^— 
or 



1 9 dv 
dr 



-u 



u\ + u\\ df B 



1 ( df B 
r V ou 2 



-e x/2 u 3 cotf 



df B df B 

U2-^ — 

ou 3 ou 2 



du 3 J 
= 0. 



(29) 



For a spherically symmetric system the coefficient of 
cot 9 in Eq. (|29|) must be zero. This implies that the 
distribution function f B is only a function of r, u\ and 
u'2 + u 2 . Multiplying Eq. (|2T)|) by mu r du, integrating over 
the velocity space, and by assuming that f B vanishes suf- 
ficiently rapidly as the velocities tend to ±00, we obtain 



rl [p{ul)]+\ P [{ul) + {A)]r d £ 



dr Lr "~ in ' 2' 
-p[{ul) + {ul)-2{ul)] 



0. 



(30) 



Multiplying Eq. (|30| by Anr 2 , and integrating over the 
cluster gives 



4np [(u 2 ) + (u 2 2 ) + (u 2 )] r 2 dr 



dv 



[(u 2 ) + (u 2 )] —dr = 0. (31) 



B. Geometric quantities 

In the following we introduce some approximations 
that apply to test particles in stable circular motion 
around galaxies, as well as to the galactic clusters. First, 
we assume that v and A are slowly varying functions of 
the radial coordinate (i.e. v 1 and A' are small). There- 
fore in Eq. (|24|) we can neglect all the quadratic terms. 
Secondly, we assume that the motion of the galaxies is 
non-relativistic, so that they have velocities much smaller 
than the velocity of the light, i.e., (u 2 ) ps (u 2 ) ph (u 2 ) <C 



and 



1. Thus, Eqs. (J24]) and ([31]) can be written as 



1 d 



2r 2 dr 



dr ) 



1 



AirGp + 4nGp 



dv 



(eff) 



2K - - / 4irr 6 p—dr = 



dr 



(32) 



(33) 



respectively, where 
r-rt 



K 



2tt P [(u\) + (u 2 ) + (u 2 )] r 2 dr, (34) 



is the total kinetic energy of the galaxies. The total 
baryonic mass of the system M B is defined as M B = 
/ dM(r) = J c I Anpr 2 dr. The main contribution to M B 
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is due to the baryonic mass of the intra-cluster gas and of 
the stars, but other particles, such as massive neutrinos, 
may also give a significant contribution to Mb- 

Now, multiplying Eq. (|32[) by r 2 and integrating from 
to r we obtain 

GMsirj^lr^-GM^ir), (35) 

where we have denoted 

M^ ff) (r) = 4tt [ p { ; ff \r')r"dr>. (36) 
Jo 

Since in hybrid metric-Palatini gravity, the quantity 
has essentially a geometric origin, we tentatively 
denote it as the geometric mass of the cluster. In the 
following we define the gravitational potential energies 
of the cluster as 



GM B {r) 



dM B (r), (37) 



f R GM^ ff \r) 

n (eff) = / i \l dMB{rl (38) 

JO r 

respectively, where R is the cluster radius. By multiply- 
ing Eq. (|35|) with rfMsfr), following an integration from 
to R, we obtain the relation 



1 

2 Jo 



R 



dv 



<>,,. ='.>:: /J "' -t / 4irr 3 p^-dr. (39) 



C. Generalized virial theorem 

Finally, with the use of Eq. (j3"3"j) , we arrive at the gen- 
eralization of the virial theorem in hybrid metric-Palatini 
gravity, which takes the familiar form 



2K + n = 0. , 



(40) 



where the total gravitational potential energy of the sys- 
tem, fl defined as 



n = n B - fi 



(e//) 



(41) 



contains a contribution term consisting of a geometric 
origin, il\ ■ 

The generalized virial theorem, given by Eq. (|40p . can 
be represented in a more transparent physical form if we 
introduce the radii Ry and R^, defined by 



Rv = Mi 



M B (r) 



dM B {r), 



(42) 



and 



R 



(eff) 



M (e//)" 


2 

/ 


'L 









R M 



{eff) 



(r) 



dM B (r), (43) 



respectively. In analogy to the geometric mass consid- 
ered above, the quantity R^ may be denoted as the ge- 
ometric radius of the cluster of galaxies in the hybrid 
metric-Palatini gravitational models. Thus, the baryonic 
potential energy VL B and the effective scalar field poten- 
tial energy fl^^ are finally given by 



n R = 



(eff) 



GM% 
Rv 



G 



M. 



(eff) 



R 



(eff) 



(44) 
(45) 



respectively. 

We define the virial mass My of the cluster of galaxies 

as 



IK 



GM B M V 
Rv 



(46) 



After substitution into the virial theorem, given by 
Eq. (|4U|) . we obtain the following relation between the 
virial and the baryonic mass of the galaxy cluster 



My 

M B 



M 



(eff) 



M%R 



(eff) 



(47) 



If My/M B > 3, a condition which holds for most of 
the observed galactic clusters, then Eq. (j47j) provides the 
virial mass in hybrid metric-Palatini gravity, which can 
be approximated as 



71/ 



(eff) 



Ry 



M B r 



(eff) ' 



(48) 



From the point of view of the astrophysical observa- 
tions the virial mass My is determined from the study of 
the velocity dispersion <r 2 of the stars and of the galaxies 
in the clusters. According to the virial theorem in hy- 
brid metric-Palatini gravity, most of the mass in a clus- 
ter with mass M to t is in the form of the geometric mass 



M 



(eff) 



, so that M 



(eff) 



M t ot- An observational pos- 



sibility of detecting the presence of the geometric mass 
and of the astrophysical effects of hybrid metric-Palatini 
gravity is through gravitational lensing, which can pro- 
vide direct evidence of the geometric mass distribution 
and of the gravitational effects associated to the pres- 
ence of the scalar field even at distances extending far 
beyond of the virial radius of the galaxy cluster. 



V. ASTROPHYSICAL APPLICATIONS 

Once the integrated mass as a function of the radius 
is determined for galaxy clusters, a physically meaning- 
ful fiducial radius for the mass measurement has to be 
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defined. The radii commonly used are either r 2 oo or 
r 50 o . These radii lie within the radii of the mean gravita- 
tional mass density of the matter (ptot) — 200p c or 500p c , 
with p c given by p c (z) — h 2 (z)3HQ/8irG, where h(z) is 
the Hubble parameter normalized to its local value, i.e., 
h 2 (z) — Vt m (1 + z) 3 + Ha, f^m is the mass density pa- 
rameter, and JIa is the dark energy density parameter, 
respectively [l8j . A pragmatic approach to the virial 
mass is to use r 2 oo as the outer boundary of the galaxy 
cluster [19]]. The numerical values of the radius r 200 are 
in the range r 200 = 0.85 Mpc (for the cluster NGC 4636) 
and r 200 = 4.49 Mpc (for the cluster A2163), so that 
a typical value for r 2 QQ is approximately 2 Mpc. The 
masses corresponding to r 2 oo and rsoo are denoted by 
M200 and Afsoo, respectively, and it is usually assumed 
that M v = M 200 and R v = r 200 'M- 



A. The weak field approximation of hybrid 
metric-Palatini gravity 

In the limit of small static gravitational fields the met- 
ric tensor can be approximated as g M „ ~ rj^ + h^, where 
t]^ is the Minkowski metric, and h^ w <C 1. The local per- 
turbation of the scalar field is denoted by ip, and therefore 
<j) k, (pQ + tp, where <po is the asymptotic value of the scalar 
field. In this approximation the evolution of the scalar 
field is described by the equation [l4j 



1 



,2 90 



Pb, 



(49) 



where ps is the mass density of the baryonic mat- 
ter, k 2 = c 3 /8t:G, and — 1/ (Gm^/c 2 ), with 
the effective mass of the field given by rrii — 
1/ [2V{4>) - V>(<f>) - 9(1 + <f>)V"( 9 )/%=t ■ 
The solution of Eq. (|4"9"j) is given by 



ip(r) 



2 GMb -2r/r v 

3 err 



(50) 



where Mb = 47r L psr 2 dr is the total baryonic mass. 
The scalar field potential can be written as V(4>) = 
V(cj> + p)^V (0 O ) + V (M p. 

In order to estimate the astrophysical effects of the 
scalar field we have to find first the explicit form of 
the effective energy associated to the scalar field p^^\ 
given by Eq. (j!?5j) . In the static case one can neglect 
all the derivatives with respect to the time, and take 



Al P = <p/rl + 
tain 



[K 2 <f>o/3\ pg- Therefore for 



we ob- 



v^ + v'iMp + ^ + ^PB 



V((j)o) + ^0 



2^ GM B 



2r/r v 



K 2 (p 



-PB, 

(51) 



where $0 = V (<po) + 1/r 2 . In the following we will as- 
sume that the baryonic matter density inside the cluster 
has a r~ 2 dependence on the distance from the cluster 
center, that is, pb oc r~ 2 , so that ps = Pbqt~ 2 , with pso 
a constant. 



B. Geometric mass and geometric radius from 
galactic cluster observations 

In the clusters of galaxies most of the baryonic mass is 
in the form of the intra-cluster gas. The gas mass density 
p g can be fitted with the observational data by using the 
following simple expression 



Pg( r ) = Po 1 + ~2 



2\ -3/3/2 



(52) 



where r c is the core radius, and po and /? are (cluster- 
dependent) constants. 

One can also assume that the pressure P g of the 
gas satisfies the ideal gas equation of state P g = 
{ksTg/ pm p )p g , where ks is Boltzmann's constant, T g is 
the gas temperature, p rs 0.61 is the mean atomic weight 
of the particles in the cluster gas, and m p is the proton 
mass [19j. This assumption is well justified physically, 
since as shown by X-ray observations, the hot, ionized 
intra-cluster gas is in isothermal equilibrium. Therefore, 
with the use of the Jeans equation [l| the total mass dis- 
tribution of the matter in the cluster can be obtained as 
a function of the gas density as [l(| E^l 



M tot (r) 



k B T, 



9 r 2 



prripG dr 



In 



Pg- 



(53) 



Taking into account the density profile of the gas given 
by Eq. (I5"2"j) . for the total mass profile inside the cluster 
we obtain the relation [l9j 



M tot (r) 



3k B (3T g 



prripG r 2 + r 2 



(54) 



According to the hybrid metric-Palatini gravitational 
model, the total mass M to t of the cluster consists of the 
sum of the baryonic mass (mainly the intra-cluster gas), 
and of the geometric mass, so that 



Mtotif) = 4tt 



Pg + P), 



1 dr. 



(55) 



Therefore Af tot (r) satisfies the mass continuity equation 
given by 



-rM tot (r) 
dr 



kKT A p a [f] + 4vrr 2 ^ e//) (r) . 



(56) 



With the use of Eqs. (|5"2"]) and (fM]) we obtain the ex- 
pression of the geometric density term inside the cluster 
as 



47T/5 



(eff) 



(r) 



3k B pT g (r 



3r 2 



AtvGpq 



pm p ( 



"c + r2 ) 2 



(l + r 2 /r 2 f P/2 ' 
(57) 
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In the cluster region where r 3> r c we obtain for p^ 
the simple relation 



(eff) 



Wl }S) (r) = 



p,m p 



1. (58) 



By power expanding the exponential in Eq. (|51[) we 
obtain 



W> e//) (r) 



c 
G 



AGM B 



3 c 2 ?v 



2^ GM B 
3 c r 



K 00 PBO 



(59) 



In the following we neglect the constant terms in the ef- 
fective geometric density, as corresponding to the cos- 
mological background. Therefore the comparison of 
Eqs. ((55} and fl5§|) gives 



*4 



k B T g 



8-kpbo l^m p 



$n = 



6nG p r c 
c 2 M B 



(60) 



Thus, from the last condition, i.e., 0o < 0, the scalar 
field potential must satisfy the condition 



V (<p ) + 1/r 2 < 0. 



(61) 



A possible potential satisfying this condition is the 
exponential potential V(cj>) = Vb exp(— 2a<p), with 
Vo and a positive constants, which gives V {4>o) = 
— 2aexp (— 2ct(f>o) < 0. Therefore for this case r v must 
satisfy the condition r v > (l/2a) exp (2a0o)- 

The geometric mass can be obtained, by using the ob- 
servational data, as 



GM^ fS \r) = 4tt 



r 2 p { f f) (r) dr 



3k B f3T g 
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inGpQ 



r dr 



pm p l + r 2 /r 2 "" ru J ( 1 + r 2/ r 2fP/ 2 
and in the limit r r c it can be approximated as 



GM, 



(eff) 



(r) 



3k B pT g AirGpor, 



pnip 



3/3„2-3/3 



3(1-/?) 



(62) 



(63) 



By using the approximation given by Eq. (1591) for the 
geometric density, we obtain 



M 



(eff) 



(r) 



c 
G 



2 GM B r 2 
3*° — J 



K 2 (p 



-pBor 



(64) 



Since astrophysical observations show that the gas rep- 
resents only a small fraction of the total mass [18l - l2ll |. 
the contribution of the gas density and mass to the geo- 
metric density and geometric mass, respectively, can be 
neglected. Therefore, for the geometric density and mass 
we obtain 



AwGp] 



(eff). 



{3k B pT, 



9 \ r -2 



(65) 



and 



GM. 



Aeff) 



(r) 



V P m P 



(66) 



(68) 



respectively. 

The form of the metric tensor component e u inside 
the cluster can be found by using the obtained mass. 
From Eq. (|5"!?|) it follows that v' satisfies the equation 
r 2 v' = 2M tot (r) ps 2M$ (r), which gives 

e" « C v r 2s , (67) 
where C v is an arbitrary integration constant, and 

3k B f3T g 

s = -. 

pGm p 

In the first approximation we can assume that the met- 
ric coefficient e _A is given by its standard general rel- 
ativistic form, e~ A « 1 — 2GM tot /r, leading to e~ A w 
1 — 6kBf3T g / pm p = 1 — 2Gs. However, this approxima- 
tion may not be valid inside the clusters of galaxies. 

Hence, in hybrid metric-Palatini gravity the metric in- 
side a galactic cluster can be directly obtained from as- 
trophysical observations. 

One may also estimate an upper bound for the cutoff 
of the geometric radius and mass, respectively. As a first 
estimate for the radius of the scalar field distribution one 
can consider the radius R<p at which the geometric den- 
sity, given by Eq. (|5"5j) . vanishes. This condition gives the 
size of the cluster as 

ksT g 



Rrh = 



(69) 



4irGpor c pm p 

The estimation of the numerical value of this radius 
requires the knowledge of po and r c , which are diffi- 
cult to be obtained from observations. An alternative 
way of defining the cutoff radius of the cluster is by 
considering the point at which the decaying geomet- 
ric density of the cluster becomes smaller than the av- 
erage energy density of the Universe. We denote the 
value of the coordinate radius at the point where the 



(cr) 



Then for r 



R 



(cr) 



two densities are equal by R 

we have (fli ) = p m iv, where p unm is the mean 

energy density of the Universe. By assuming p un t v — 
p c = 3H 2 /8ttG = 4.6975 x lCT 30 ^ g/cm~ 3 , where 
H = 50h 50 km/Mpc/s [H, 

we obtain 



R 



(cr) __ 

I 3ks/3T g ^ 



\pm p Gp c 



5keV/ 



Mpc. 



(70) 



The total geometric mass = yR^j extend- 



r>(cr) . 

mg up to R\ is 



r> = 4.83 x lO-W^f \- 



M 



V 5keV J 



Mr. 



(71) 
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This value of is consistent with the observations 

of the mass distribution in the clusters of galaxies fl9j . 
However, according to the hybrid metric-Palatini gravi- 
tational model, we predict that due to the presence of the 
effective scalar field the geometric mass and its effects do 
extend beyond the virial radius of the clusters, which is 
of the order of only a few Mpc. 



C. Radial velocity dispersion in galactic clusters 

In terms of the characteristic velocity dispersion a\ the 
virial mass can also be expressed as [21| 



M v = -<j{Rv. 



(72) 



Taking into account that the velocity distribution in the 
cluster is isotropic, we have (u 2 ) = (u 2 ) + (u 2 ,) + (u 2 ) = 



3(u 2 ) = Z<J 2 , where a 2 is the radial velocity dispersion. 
o i and ay are related by ia\ = a 2 .. 

The radial velocity dispersion relation for clusters of 
galaxies in hybrid metric-Palatini gravity can be derived 
from Eq. (|3U|) . Since the velocity distribution is isotropic, 
we obtain first 



dr 



\pa r 



1 dv 
2'dr =0 - 



(73) 



Inside the cluster in the first order of approximation 
the condition e~ A « 1 holds. In the limit of small ve- 
locities the modified field equation, Eq. (IM1) can be inte- 
grated to yield 



rV 



2GM <t> (r) + 2GM B (r) + 2C, 



(74) 



where C is an arbitrary constant of integration. Since 
from Eq. (|73|) we have v' = —(2/p)d (per 2 ) /dr, it follows 
that in hybrid metric-Palatini gravity the radial velocity 
dispersion of the galactic clusters satisfies the differential 
equation 



(75) 



with the general solution given by 



<{r) = — 



GM^r 1 ) 



P 

GM(r' 



■p(r') 



p(r') 
3^ 



dr 1 



, (76) 



where Ci is an integration constant. 

In the following we consider a simple case in which the 
density p of the normal matter inside the cluster has a 
power law distribution, given by 



p(r) = p r~ 



(77) 



with pq and 7 ^ 1,3 positive constants. The cor- 
responding baryonic matter mass profile is Mb(t) — 



inpor 3 7 / (3 — 7). The geometric mass, given by 
GM^ fa (3k B f}Tg/ptm p )r = q r, where q Q = 
(3fc_B (3T g I prrip) , is linearly proportional to r, as has been 
shown in the previous Section. 

Therefore, for 7 7^ 1,3, we obtain the following expres- 
sion for the velocity dispersion, 



2ttG Po 



,2- 7 , C 1 . Ck 



-r 2 ^ + -- +—r\ (78) 



7 (7 - 1) (3 - 7) 7 + 1 r po 

For 7 = 1, we find 

C C 
<?r (r) = 90 + 7T - 2nGp r In r + — r, (79) 
2r p 

while for 7 = 3 we obtain 

— A-.I)4 + £I + ^. ,80, 
4 / r 4 4 r po 



a r ( r ) = y - nG Po ( In r 



For clusters of galaxies the observed data for the ve- 
locity dispersion are usually analyzed by assuming for 
the radial velocity dispersion the simple form <r 2 (r) = 
B/(r + b), with B and b constants. For the density of the 
galaxies in the clusters the profile p(r) = A/r (r + a)~ , 
with A and a constants, is used. The observational data 
are then fitted with these functions by using a non-linear 
fitting procedure plj . For r -C a, p(r) s» A/r, while for 
r S> a, p(r) behaves like p(r) ps A/r 3 . Therefore the 
comparison of the observed velocity dispersion profiles 
of the galaxy clusters and the velocity dispersion pro- 
files predicted by hybrid metric-Palatini gravity provides 
a powerful method for observationally discriminating be- 
tween the different modified gravity theoretical models. 



VI. DISCUSSIONS AND FINAL REMARKS 

Modern astrophysics and cosmology have to face two 
major difficulties: the dark energy and the dark mat- 
ter problems, respectively. One promising approach to 
improve our understanding of these issues is to modify 
gravity at large galactic and cosmological scales. In par- 
ticular, the hybrid metric-Palatini gravitational theory 
can challenge the need for dark matter and dark energy. 
In the framework of this theory, cosmological models that 
account for the late time acceleration of the universe do 
exist [l5j], as well as viable models that pass all the re- 
quired Solar System tests [14 1. 

Moreover, as shown in the present paper, in the frame- 
work of hybrid metric-Palatini-type modified theories of 
gravity, the galactic dynamics of massive test particles 
may also be understood without the need for dark mat- 
ter. We have analyzed the "dark matter" problem by 
considering a generalized version of the virial theorem. 
The virial theorem was obtained by using a method based 
on the collisionless Boltzmann equation. The additional 
scalar field terms present in the modified gravitational 
field equations give an effective contribution to the grav- 
itational energy, which at the galactic/extra-galactic level 
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acts as an effective mass, piaying the roie of the "dark 
matter". The total virial mass of the galactic clusters 
is mainly determined by the effective mass associated to 
the new scalar field terms, which can be called geomet- 
rical mass, since its intrinsic origin is geometrical. The 
presence of this term may explain the well-known virial 
theorem mass discrepancy in clusters of galaxies [l| . 

In the framework of the hybrid metric-Palatini gravita- 
tional theory we have also shown the existence of a strict 
proportionality between the virial mass of the cluster and 
its baryonic mass, a relation which can also be tested ob- 
servationally. Since galaxy clusters are "dark matter" 
dominated objects, the main contribution to their mass 
comes from the geometric mass mW^\ so that with a 

very good approximation we have My M to t- 

Therefore the virial theorem gives immediately the fol- 
lowing mass scaling relation 

R { I ff) 

My « Mb — 5 — • (81) 
riy 

This equation shows that in the framework of hybrid 
metric-Palatini /(X)-type modified gravity the virial 
mass is proportional to the baryonic (normal) mass of 
the cluster, and that the ratio of the total mass and of 
baryonic mass is determined by a purely geometric quan- 
tity, the geometric radius R$. Hence the geometric ra- 
dius of the cluster can be determined from observations, 
once the virial and baryonic masses and the virial radius, 
respectively, are known. 

By assuming that the cluster extends up to a value of 
the radius given by R^ rs Rl ° , we obtain the following 
relation between the virial and the baryonic mass of the 
cluster 

For a cluster with gas temperature T g = 5 x 10 7 K, f3 = 
1/2 and Ry = 2 Mpc we obtain My ~ 32Mb, a relation 
which is consistent with the astronomical observations 
19]. 

One of the important, and observationally testable, 
predictions of the hybrid metric-Palatini gravitational 



"dark matter" model is that the geometric masses asso- 
ciated to the clusters, as well as its gravitational effects, 
extend beyond the virial radii of the clusters. Observa- 
tionally, the virial mass My is obtained from the study of 
the velocity dispersions of the stars in the cluster. Due 
to the observational uncertainties, this method cannot 
give a reliable estimation of the numerical value of the 
total mass M g + M^^ in the cluster. However, a much 
more powerful method for the determination of the total 
mass distribution in clusters is the gravitational lensing 
of light, which may provide direct evidence for the grav- 
itational effects at large distances from the cluster, and 
for the existence of the geometric mass. The presence of 
hybrid metric-Palatini modified gravity effects at large 
distances from the cluster, and especially the large ex- 
tension of the geometric mass, may lead to significantly 
different lensing observational signatures, as compared 
to the standard relativistic/dark matter model case. The 
bending angle in the hybrid metric-Palatini gravity mod- 
els could be larger than the one predicted by the stan- 
dard dark matter models. Therefore, the detailed obser- 
vational study of the gravitational lensing could discrim- 
inate between the different theoretical models introduced 
to explain the motion of galaxies ( "particles" ) in the clus- 
ters of galaxies, and the standard dark matter models. 

To conclude, the generalized virial theorem in hybrid 
metric-Palatini gravity might be an efficient tool in obser- 
vationally testing the viability of this class of generalized 
gravity models. 
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